1 



o 
o 
o 

(N 

< 

(N 



Determinations of \Vub\ with inclusive techniques at LEP 

M. Battaglia''* 

''Dept. of Physics, University of Helsinki, 
FIN-00014 Helsinki, Finland 

The charmless semileptonic decay B branching fraction has been measured, using inclusive techniques, by the 
Aleph, Delphi and L3 experiments at LEP. The average of their results is 

BR(6 ^ = (1.74 ± 0.37 (stat.+exp.) ± 0.38 (6 ^ c) ± 0.21 [b ^ u)) x 10"^ 

From this result the value of the \ Vub \ element in the CKM mixing matrix has been derived, using OPE predictions, 
obtaining: 

\Vub\ = (4.13lJ5:^? (stat. +det.)+°;*| (b ^ c syst.)+°;^* (b ^ u sys.) ±0.02 (rs) ± 0.20 (HQE)) x 10"^ 
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1. INTRODUCTION 

The accuracy in the determination of the \Vub\ 
element in the Cabibbo-Kobayashi-Maskawa mix- 
ing matrix 0] plays an important role in the study 
of the unitarity triangle and in the related tests of 
the Standard Model (SM). In particular, a non- 
vanishing \Vub\ is essential to preserve the possi- 
bility to describe CP violation within the SM and 
its value places a direct constraint on the magni- 
tude of the CP violating phase /?. 

The first determination of the magnitude of 
\Vub\ was obtained from the yield of leptons 
produced with momentum above the kinematic 
limit for b Xciv transitions, first reported by 
Cleo [H and soon confirmed by Argus Q . How- 
ever this method is sensitive to only ~ 10% of the 
inclusive charmless semileptonic (s.l.) yield, and 
the extraction of \Vub\ is subject to a large model 
dependence. 

More recently, exclusive B -niv and B — > piv 
decays have been measured by Cleo The 
determination of IKitj from exclusive s.l. decays 
has still a significant model dependence. First lat- 
tice estimates of the relevant form factors indicate 
that a significant reduction of these uncertainties 
may be expected in the future. 
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The extraction of |V^b| from the distribution 
of the invariant mass Mx of the hadronic sys- 
tem recoiling against the lepton pair peaked, for 
h — > Xy_lv^ at a significantly lower value than for 
h X(.lv was proposed several years ago [|| , and 
it has recently been the subject of new theoretical 
calculations . If 6 — > m transitions can be dis- 
criminated from the dominant 6 — *■ c background 
up to Mx — M{D), this method is sensitive to 
~ 80% of the charmless s.l. B decay rate. Fur- 
ther, if no preferential weight is given to low mass 
states in the event selection, the non-perturbative 
effects are expected to be small and the OPE de- 
scription of the transition has been shown to be 
accurate away from the resonance region. 

The experimental challenge comes from the re- 
quirement to isolate the 6 — > w contribution to 
the s.l. yield from the ~ 60 times larger 6 — > c 
one while ensuring a uniform sampling of the de- 
cay phase space to avoid biases towards a few 
exclusive low-mass, low-multiplicity states, such 
as tt£9 and piD. In principle, this method is 
well suited for the LEP experiments, where the 
recorded statistics is not sufficient for studying 
the exclusive decay modes with good accuracy 
while the significant boost of the B hadrons, the 
separation of b and b decay products in opposite 
hemispheres, and the good secondary vertex re- 



2 



construction capabilities make the study of inclu- 
sive decay B decays possible. The analysis tech- 
niques adopted by the Aleph |^ , Delphi j|] and 
L3 Collaborations are based on the obser- 
vation that b XuiP decays can be inclusively 
discriminated from b Xciv by exploiting the 
differences in the invariant mass and kaon con- 
tent of the secondary hadronic system, in the de- 
cay multiplicity and in the decay vertex topology. 
These features have been used differently in the 
three analyses, resulting in determinations of the 
charmless s.l. branching fraction obtained from 
samples with varying efficiency and purity and 
with systematic uncertainties that are only par- 
tially correlated. 

2. EXPERIMENTAL RESULTS 

The analyses consist of three main steps: i) in- 
clusive selection of s.l. B decay candidates, ii) 
definition of a subsample enriched in 6 — > w tran- 
sitions and iii) measurement of the charmless s.l. 
branching fraction. All experiments started from 
a sample of events containing an identified elec- 
tron or muon selected from a combined data set of 
~ 8.2 M hadronic Z° decays. Aleph and Delphi 
also imposed 6-tag criteria to reject light quark 
and charm backgrounds. Aleph used two neu- 
ral networks (NN) to separate charged and neu- 
tral B decay products from fragmentation par- 
ticles, reconstructed the secondary hadronic sys- 
tem and the B rest frame and computed twenty 
b ^ u discriminating kinematical variables in this 
frame. Delphi also reconstructed the secondary 
hadronic system mass and the B rest frame, but 
using a particle likelihood variable and an itera- 
tive topological reconstruction procedure and es- 
timated the lepton energy in the B rest frame 
i?| for each event. L3 adopted a consecutive 
cut analysis based on the kinematics of the two 
most energetic hadrons in the same hemisphere 
as the tagged lepton. All the three experiments 
observed a significant excess of events with the 
characteristics expected for b — > X^iv decays. 

Aleph combined the selected discriminating 
variables by means of another NN to obtain a 
global discriminant variable NNf,„ (see Figure |^) . 
The number of & — > X^iv candidates in the data 
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Figure 1. The b u discriminating NN output 
for the Aleph analysis, a) and c) include only the 
simulated backgrounds while b) and d) give the 
comparison with the data after including 6 — > m 
s.l. decays corresponding to the fitted branching 
fraction 

was extracted by a binned likelihood fit to the 
NNfeu and converted into the charmless s.l. BR. 
Delphi divided the reconstructed events in four 
classes on the basis of the Mx value and of en- 
richment criteria based on the relative position of 
the lepton w.r.t. the secondary vertex and the 
presence of tagged kaons. Of these classes the 
Mx < 1.6 GeV/c^ - b ^ u enriched one was ex- 
pected to contain almost 70% of the s.l. b ^ u 
decays while the other classes, depleted of sig- 
nal events, were used to monitor the background 
modeling in the simulation (see Figure ^). The 
ratio |K,f)|/|V^cfc| was extracted, together with the 
overall data/MC normalization, from the fraction 
of b ^ X^lv candidates observed in the data us- 
ing a two parameter likelihood fit to the number 
of events in each of the four classes and to their 
E"^ distributions. Finally L3 extracted the charm- 
less s.l. B branching fraction by counting the ex- 
cess of events over the estimated background after 
having applied their selection criteria. 

Starting from a natural signal-to-background 
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Figure 2. Background subtracted E'^^ distribu- 
tions for the Delphi analysis: the b u enriched 
decays with Mx < 1.6 GeV/c^ (upper plot) and 
h ^ u depleted decays with Mx < 1-6 GeV/c^ 
(lower plot) . The shaded histograms show the ex- 
pected distribution for signal b ^ u s.l. decays 
normalized to the amount of signal corresponding 
to the fitted iKibl/lV'cbl value. 
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Figure 3. Output of a NN combining the discrim- 
inating variables of the L3 analysis. The amount 
of signal in the simulation has been scaled to the 
measured BR. This NN has only been used as 
consistency check. 



ratio, S/B, of about 0.02, Aleph obtained S/B 
= 0.07 with an efficiency e = 11%, Delphi S/B 
= 0.10 with e = 6.5% and L3 S/B = 0.16 with 
e — 1.5%. The results are summarized in Ta- 
ble ^ Several consistency checks were performed 
by the three experiments. Aleph explicitly used 
vertexing variables, Delphi performed a search 
for fully reconstructed B — > tt£D, plD decays and 
L3 validated the excess of events by constructing 
a 6 ^ u discriminating NN (see Figure ||) . 

3. AVERAGE BR(6 ^ XJv) AND 
EXTRACTION OF |K&| 

The three measurements of BR(6 X^^v) 
have been averaged using the Best Linear Un- 
biased Estimate (B.L.U.E.) technique [Hj. This 
technique provides with an unbiased estimate 
BRlep that is a linear combination of the differ- 
ent measurements BRj corresponding to the min- 
imum possible uncertainty a: 



BRlep — 



EliE,tiBRi(ii;'^). 



(1) 



with (j^ = 



.^^3 — where E is the error 

matrix including the off-diagonal terms giving the 
correlations between pairs of measurements. 

The sources of correlated systematics belong 
to both the description of background b c 
and to the modeling of signal b ^ u transitions. 
The differences in the analysis techniques adopted 
by the three experiments are reflected by differ- 
ences in the sizes of the systematic uncertain- 
ties estimated from each common source. Impor- 
tant common systematics are due to the charm 
topological branching ratios and to the rate of 
D decays. Aleph and L3 are also sensi- 

tive to the uncertainties in the b fragmentation 
function due to the use of kinematical variables 
for enriching in 6 ^ X^tv. The Delphi result 
depends on the assumed composition in b hadron 
species due to the use of kaon anti-tagging to re- 
ject b ^ c, thus rejecting also Bg and At de- 
cays. The signal b —>■ u systematics have been 
grouped in inclusive model and exclusive model 
and assumed to be fully correlated. The first cor- 
responds to the uncertainty in modeling the kine- 
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Table 1 



The results on BR(6 Xu^v) from the LEP experiments with the sources of the statistical, experimental, 
model uncorrelated and model correlated uncertainties. 



Experiment 


BR(6 ^ XJv) 


(stat.) 


(exp.) 


(uncorrelated) 


(correlated) 


Aleph §1 
Delphi § 
L3 @ 


(1.73 
(1.69 
(3.30 


± 0.48 
± 0.37 
± 1.00 


± 0.29 
± 0.39 
± 0.80 


± 0.29 IZl) 
± 0.18 {f,'^ IZl) 
± 0.68 (±"l68 b^cj 


±0.47 (±°:« izD) xio-3 

±0.42 {f,%l IZD) xlO-3 

±1-40 {f,f, IZD) xio^^ 



matics of the 6-quark in the heavy hadron. It 
has been estimated from the spread of the results 
obtained with the ACCMM model |0, a shape 
function, describing the distribution of the light- 
cone residual momentum of the heavy quark in- 
side the hadron [|l^-|l^ and the parton model ||l^ 
in the Aleph and Delphi analyses and from the 
uncertainties in the single tt and the lepton energy 
spectra for L3. The exclusive model uncertainty 
arises from the modeling of the hadronic final 
state in the b — > Xuiv decay. These uncertain- 
ties have been estimated by replacing the parton 
shower fragmentation model in JETSET |l^ with 
the fully exclusive ISGW2 ||l|| model by Aleph 
and Delphi and by propagating a 100% uncer- 
tainty on the B tt£D rate by L3. Using the 
inputs from Table |l|, the LEP average value for 
BR{b ^ XJD) was found to be: BR(6 XJD) 
= (1.74 ± 0.37 (stat.-Hexp.) ± 0.38 (6 ^ c) ± 0.21 
(6 u)) X 10-^ = (1.74 ± 0.57) x 10"^ ^j^h a 
confidence level for the combination of 0.723 [l9| . 

The value of the \Vub\ element has been ex- 
tracted by using the following relationship de- 
rived in the context of Heavy Quark Expan- 
sion MM 



I Kb I - .00445 ( 



BR(6 Xulv) 1.55ps 



0.002 



n 



)2 X A (2) 



with .4 = (1 ± 0.020(QCD) ± 0.035(mb)) where 
the value to;, = (4.58 ± 0.06) GeV/c^ has been 
assumed [p2| . The uncertainties have been con- 
voluted together assuming them to be Gaus- 
sian in BR(fe — s- Xulv)^ with the exception of 
the theoretical uncertainty on A assumed to 



= 4.13 
+ 0.42 
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(stat+exp.) 

+ 0.24/h 
- 0.25 

± 0.20 (HQE) 
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Figure 4. The probability density function for 
I Kb I corresponding to the LEP average value of 
BR(fe Xu^v) with the intervals corresponding 
to the 68% C.L. and 95 % C.L. indicated. 



be Gaussian in \Vub\- The resulting probability 
density distribution is shown in Figure ^ and 
gives iKbl = (4.13t°:?|) X 10-3 at 68% C.L. 
and \Vub\ = {'^■Vit\jl) X 10-3 at 95% C.L. The 
part of this function in the negative, unphysical 
region is negligible, corresponding to only 0.12%. 
By repeating this procedure separately for each 
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systematic, the detailed result for the 68% C.L. is: 
\Vub\ = (4.13t°:^?(stat. + det.)+o:48(b ^ c syst.) 
to.25(b ^ u sys.) ±0.02(Tf,)±0.20(HQE)) x IQ-^. 



4. DISCUSSION AND CONCLUSIONS 

The LEP analyses have demonstrated the feasi- 
bility of an inclusive determination of the charm- 
less s.l. B branching fraction by discriminating 
h ^ u from 6 — > c decays on the basis of the 
mass, multiplicity and kaon content of their sec- 
ondary hadronic system and of their decay topol- 
ogy. Differently exploited by three of the LEP ex- 
periments, these event characteristics have been 
used to obtain clear signals for the decay, to mea- 
sure its branching fraction and to derive a LEP 
combined value of = (4.13tEJ;^^) x 10"^ at 
68% C.L. with the relative uncertainty due to the 
6 — > u model below 10%. This result agrees with 
the recent Cleo determination Q using the ex- 
clusive B — > piD decay branching fraction giving 



1Kb I = (3.25l^;^4) X 10^3 where the uncertainty 
is dominated by a 17% model systematic mostly 
uncorrelated with that of the LEP measurement. 
The agreement between the inclusive and exclu- 
sive determinations, as in the case of \Vcb\^ is en- 
couraging as a test of the underlying theory as- 
sumptions and to control possible violations of 
quark-hadron duality in semileptonic B decays. 
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